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Abstract: We explore in the supergravity context the possibility that a Higgs scalar 
may drive inflation via a non-minimal coupling to gravity characterised by a large 
dimensionless coupling constant. We find that this scenario is not compatible with the 
MSSM, but that adding a singlet field (NMSSM, or a variant thereof) can very naturally 
give rise to slow-roll inflation. The inflaton is necessarily contained in the doublet Higgs 
sector and occurs in the D-flat direction of the two Higgs doublets. 
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1. Introduction 

The idea that a large non-minimal coupling of scalar fields to gravity might play an 
important role in inflation is clearly worth exploring. (For an early example see Ref. |l|]) . 
In particular, there has been recent interest in the economical possibility that 

the standard model (SM) Higgs H can be the relevant scalar; the action 5* is 



where R is the Ricci scalar, and £ is a dimensionless coupling constant, assumed pos- 
itive, 1 and Csm is the Lagrangian of the SM, generalized to a nontrivial background 
spacetime obtained by replacing the flat-space metric r]^ v with the curved space metric 
g^y. Classically, for large enough values of £ and h, viz., 



(where H = (0,h) T ), it transpires that the scalar potential is nearly flat and the 
standard slow-roll approximation is possible. The issue of whether this classical analysis 

^^Hcre, the action has been expressed in the Jordan frame; one may transform to the Einstein frame 
by a conformal transformation. 




(1.1) 



til 



2 > M 2 P » h 2 , 



(1.2) 



remains valid as an effective quantum field theory for field values in the range Eq. ( |1.2|) 
is a matter of debate at present. On the one hand, there are some || [7] who argue that, 
as an effective field theory, the full action Eq. (|1.1| ) must include higher dimensional 
operators, such as £, 2 {H^ H) 6 /Mj>, so that one would expect the theory to work up to 
an energy scale A = Mp/£, which is far below A = Mp/y^, the regime suggested by 
Eq. QT3). 

On the other hand, others argue the effective field theory works within the regime 
in Eq. (|1.2j ) for a variety of reasons. One such argument fl], f|, ||, |[ is that loops 
involving virtual Higgs are suppressed because their propagators are modified from 
ordinary perturbation theory for large values of the Higgs field. The treatment of 
gravity as a classical background may also be challenged. One would normally expect 
that the feedback on the metric would become strong on scales of the order of Mp, but 
one may hope [|l], |2[ that, in some models, the the energy density associated with Higgs 
remains small compared to Mp, even for h > O(Mp), provided it is not very much 
larger, say h < O(10Mp). On the other hand, Higgs-graviton scattering would seem to 
become strong or to violate unitarity for energies E ~ Mp /£ |], . We will not address 
this issue directly in this paper, but we will return to it again in our conclusions. 

The occurrence of such non-minimal couplings of matter fields to the background 
curvature is a fact of life in the SM, because the notion of minimal coupling (£ = 0) is 
only valid classically. In quantum field theory, the renormalized £ runs, i.e, is a function 
of a scale parameter (£ = for which £ = is not a fixed point [12], |j~3[| . Thus, 



the "minimal theory" is at best an approximation to cases where the contributions of 
terms like the £ term above may be neglected. 

Whether or not perturbation theory remains viable in the scenario, the scales of 
interest are large compared to the electroweak unification scale, so that the SM becomes 
unnatural, suffering also other apparent shortcomings such as unsatisfactory grand 
unification and lack of a suitable Dark Matter candidate. The latter issue is addressed, 



for example, in Refs. |I4| [L5[ by adding a singlet to the standard model. More attractive 
alternatives to the SM above the TeV-scale are supersymmetric extensions of SM, which 
overcome the three most glaring deficiencies of the SM mentioned. An obvious question 
is whether the inflationary scenario described above survives in such supersymmetric 
models, and, if so, what their properties are. In this paper, we shall consider the MSSM 
and a simple extension with an additional singlet field (NMSSM), with which we will 
assume the reader is familiar (for a review see Ref. [0). 

An outline of the remainder of this paper is as follows: in the next section, we 
review the standard form of supergravity and then consider the addition of a non- 
minimal interaction analogous to the above. In Section 3, we apply this formalism to 
the MSSM , following in Section 4 with the NMSSM. Section 5 contains an analysis of 



-2- 



the slow roll parameters in the one successful scenario that we identified. In Section 6, 
we restate our conclusions and discuss issues for future work, commenting on some 
dramatic differences from the SM to be anticipated for radiative corrections resulting 
from supersymmetry. 



2. Non-minimal couplings in supergravity 

Although particle physics generally ignores gravity in generalizing to supersymmetry, 
its inclusion is mandatory in the present context. Thus, the natural starting point is 
not global supersymmetry but the modification of the Lagrangian for supergravity 2 
coupled to a multiplet of chiral superfields $: 



C 



6 / d 2 QS 



P-^-(£> 2 -8P)$ f $ + P($) 



h.c. + gauge terms. 



(2.1) 



where £ is the vierbein multiplet, R is the curvature multiplet containing the scalar 
curvature in its 2 component, and P($) is the superpotential. For our purposes, we 
will not need the explicit form of the gauge field contributions. 

To extend Eq. (|1.1|) to the supersymmetric case we replace Eq. (|2.1|) by 



£ x = -6 



d 2 QS 



R + X($)R - i(£> 2 - 8R)&$ + P($) 



+/i.c.+ gauge terms (2.2) 



where X($) is quadratic in $, with coefficients which are dimensionless coupling con- 
stants. We shall discuss examples below, but an immediate consequence of SUSY is 
that the gauge invariant function A($) is a function of the chiral superfields $ only and 
not $t. A second consequence is that, by the non-renormalization theorem, the various 
monomials in X(<&) will (in a similar manner to those in the superpotential P($)) 
not be generated through radiative corrections if absent from the classical action (in 
contrast to the case of the SM £ from Eq. ( |1.1|) ). 

The scalar potential in the Einstein frame for the theory defined by Eq. (|2.1|) is 
given by 



V(0,0* 



[Gi(G- v 



ij*Gj* 



3}+V D ( 



e K [{K-^.DiPiDjPy - 3PP*] + V D {(j>, , 



(2.3) 



where G = K + ln(PP*), K((/>, 0*) is the Kahler potential, and the D-terms Vd( 
will be discussed below. We will take K to have the canonical form 

K = -31n(-fi/3), 



(2.4) 



2 Our notation and conventions follow Wess and Bagger [0, and from now on we set Alp = 1. 
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where 



3. 



(2.5) 



In Eq. (g^), D { P = d { P + K { P, etc. 

One can show that the effect of the non-minimal couplings on the scalar potential, 
generalising from Eq. (|2.1|) to Eq. (p.2|), is simply to replace Q in Eq. (|2.5|) by 



3--(X(0) + / i .c). 



(2.6) 



Let us turn to the .D-term. In general, this takes the form (for a simple group with 
gauge coupling g) 



V D = 9 -Ref^ b 1 G i {T a ) ij <P 3 G k {T h ) kl< p h 



(2.7) 



where f a t, is associated with the kinetic energy of the gauge field and, therefore, must 
be a holomorphic function of (pi, and where (pi transform according to a representation 
T a of the gauge group. Noting that 



diP(T a ] 



«9,X(T a ; 







(2- 



by gauge invariance, and assuming for simplicity the canonical form f a b = 5 a b, we find 



V D 



Thus in the MSSM, the D— terms are 



2 0*' 



(2.9) 



(2.10) 



where Hi and H2 are the two Higgs doublets, and f are the Pauli matrices. 

The behavior of the potential depends in detail on the choice of superpotential 
P{(p), but we will be interested in directions in field space in which the potential is 
approximately constant for large values of the Higgs fields. 



3. MSSM 

Let us begin by considering the effect of Eq. (|2.6|) on the scalar potential for the MSSM. 
The unique possibility for X in this case is 

X = xH^, (3.1) 
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where x is constant; and, dropping terms involving fields that do not appear in X, the 
most general possibility for the superpotential P is 

P = A + iiH x H 2 , (3.2) 

where A and /i are constants. So P contains no dimensionless couplings. One cannot 
in general choose both \i and x to be real; however, we will ignore this potential source 
of CP-violation in what follows, since, as we shall soon see, this case does not yield a 
suitable inflationary regime. Without loss of generality, we may choose x > 0. 
Let us calculate the potential in the electromagnetism-preserving direction 

H 1 = (h 1 ,0) T ,H 2 = (0,h 2 ) T , (3.3) 

where we may assume h 1}2 are real. If, as usual, we define h x = h cos (3 and h 2 = h sin (3, 
we can then write 

P = A + ///i 2 sin(2/3)/2, X = X h 2 sin(2/?)/2. (3.4) 

We shall be interested in the case when x ^> 1, but h 2 <C 1, so we must have x sin(2/3) > 
to avoid a singularity in Q x . Since x > 0, we must take < (3 < 7r/2(mod ir). We 
find 

o ( 6A + sin W + 12 X^ 2 (3A + fxh 2 sin 2(3) - Yltfh 2 

{3 X h 2 sin 2(3 + 6- 2h 2 ) 2 (3 X 2 h 2 - 2 X h 2 sin 2(3 + 4) + D ^ ' ' 

where 

9(^ + ^ 2 )^cos 2 2/? 
D 2(3x/i 2 sin2/5 + 6-2/i 2 ) 2 ' 1 J 

It is easy to see that, for fixed (3 in the regime 

X h 2 > 1 > /i 2 (3.7) 

the potential is dominated by the D-term, 

V D ~ g -±l- cot 2 2/5, (3.8) 

except very near the .D-flat direction, more precisely, for sin 2(3 3> 2/(xh 2 )- It is clear, 
however, that slow-roll conditions (to be reviewed below) will not be satisfied in these 
circumstances, because 

1 dV D 



Vd dp 



(3.9) 



sin 4,3| 



cannot be made small. 
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On the other hand, in the .D-flat direction where tan/? = 1, both Vd and dVo/dfi 
vanish, and d 2 Vd/d 2 P > 0. Thus, this is a minimum stable against fluctuations in (3. 
Then, in the regime described in Eq. ( |3.7| ), we have 



V 



2[3A 2 + 2fi X h 2 {2fih 2 + 3A)] 
3?^ 



(3.10) 



or for A = 0, V — —8fi 2 /(3x 3 h 2 ). In either case, V —>■ through negative values, and 
so is unsuitable for an inflationary scenario. 

4. NMSSM 

Evidently, the failure of this MSSM to emulate the result of a nonminimal coupling in the 
nonsupersymmetric case is because of the absence of a cubic term in the superpotential, 
leading to the absence (in the .D-flat direction) of an analog to the dimensionless self- 
interaction A(0V) 2 - (Of course inclusion of the usual Yukawa cubic terms in P would 
lead to quartic terms in V, but we need quartic terms that contain only fields that 
occur in X in order to produce a suitable flat potential). It is therefore natural to turn 
to the NMSSM, which clearly can provide a cure to this behavior. Thus we consider 
the superpotential 



where S is a gauge singlet. There are then two natural choices for X (</>), to wit 




(4.1) 



X = X HiH 2 or 
X = X S 2 . 



(4.2) 
(4.3) 



Note that in the case of Eq. ( |4.2j ) we can choose both \ an d A to be real. 

If we choose Eq. ( |4.2|) and consider the case when S is small but h is large, 



S = 0,H 1 = (h 1 ,0) T ,H 2 = (0,h 2 ) T 



(4.4) 



then we find 



9(2A 2 /i 4 sin 2 2(3 + {g 2 + g' 2 )h^ cos 2 2(3) 
~ 2(3xh 2 sin 2(3 + 6 - 2h 2 ) 2 

which, taking the the Eq. ( p.7|) limit, becomes 



(4.5) 




(4.6) 
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Figure 1 
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Figure 1: Plot of V against h 

which is (of course) minimized for tan/3 = 1. For this value of tan/3, Eq. (|4.5| ) gives 

9X 2 h 4 



V 



(3 X h 2 + 6-2/i 2 ) 2 ' 



(4.7) 



a smoothly increasing function of h 2 , precisely the behavior that we can expect to yield 
inflation. This is the supersymmetric analog of the model of Ref. 0. We plot V in 
Figure 1, for £ = 5 x 10 4 and A = 1. 

On the other hand, if we choose Eq. ( |4.<j|) and h small, then we find 3 



V = 36\p 2 



[ X (S 2 + S* 2 )-2}S 2 S* 2 



[3 X (S 2 + S* 2 ) - 2SS* + 6} 2 [ X {S 2 + S* 2 ) - 6 X 2 SS* - 2] ' 

Writing S = se td we have 

y = 9 | 2| ( X s 2 COS26 - l])s* 

lP 1 (3xs 2 cos 26 - s 2 + 3) 2 (xs 2 cos 26 - 3x 2 s 2 - 1] 

which in the limit X s " 2 ^ I ^ s 2 gives 



(4.8) 



(4.9) 



V 



\P 



3x 3 cos 26 



(4.10) 



3 In principle, we cannot choose both p and \ reai - For our purposes, it is convenient to take x rea -l 
but p possibly complex. 
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Thus for fixed 6, V approaches a constant, but if xcos2# > then this constant is 
negative. If we choose x cos 2$ < in order to obtain V > then we see from Eq. ( |4.9|) 
that as s reduces, then V increases, and we eventually encounter a pole in V, when 

XS 2 cos29 = -1 + s 2 /3 (4.11) 

Thus this case is not suitable for inflation. 



5. Slow-roll analysis 

Let us pursue the more promising case of Eq. fl4.7|) in more detail and calculate the 
slow-roll parameters. Evidently these relate to fluctuations and so we must consider 
the kinetic energies of the relevant fields. It is easy to show using Eq. ( |2.6|) that in 
general the scalar kinetic energy takes the form 

£ kin = -e{K x ) ir ,d m ^d m <P) (5.1) 

where e is the vierbein and 

K x = -31n(-O x /3). (5.2) 



Let us now briefly review the slow- roll paradigm (see for example Ref. |18]). Suppose 
the action that describes the inflaton 6 takes the form 



S — id xe 



Then neglecting its spatial variation, the equation of motion for is 



(5.3) 



mi + 3 w + \rm 2 + ^ = ( 5 - 4 ) 

where 

H 2 = \ Q/(0)0 2 + V {<(>)) ■ (5.5) 



Then the slow-roll parameter e is given by 



H 1 fV'^ 2 



2H 2 2f \ V 



(5.6) 



where the latter approximation is valid if f{4>)4> an d /'(0)0 2 can both be safely neglected 
in Eq. fl5l]). 
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The number of e-folds before inflation ends is given by 



r H r fV 
N= T d^« / J —d4 (5.7) 

"''/'end <P "'^cnd 

and the second slow-roll parameter r\ is given by 

V" 

V = jy- (5-8) 

In the slow-roll regime, e, \rj\ <C 1 and slow-roll ends when e(0 en d) ~ 1. 

We now apply this formalism to the NMSSM doublet case. We first suppose that 
tan/3 = 1, that is we consider the .D-flat direction. From Eqs. ( ^6|) , (|3.3|) , (|4.2| ), ( [5.1|) 
we find that for large x ( as defined in Eq. ( |3.7|) ) the dominant contribution to the 
inflaton kinetic energy term is 

L kin = 3e(^\\ (5.9) 



Armed with this result we immediately obtain (given Eq. (|3.7| )) from Eqs. Q5.6|) , 
~>J\) and (|5~B that 

16 

(5.10) 



8 V ' 

V = ~. (5.12) 
Xh 2 



We see that the slow-roll conditions e,\r]\ <C 1 are satisfied given Eq. ( |3.7| ). Moreover 
it is easy to verify that the approximations made in neglecting f(4>)4> and /'(</>) 2 in 
Eq. ( |5.4| ) are well justified. The results are similar to those of Ref. H. For e(/i en d) ~ 1, 
and using x ~ 5 x 10 4 as in Ref. |J, we obtain h end ~ 0.7 x 10~ 2 Mp, in the neighborhood 
of the GUT scale. 

Suppose that tan /3 1 at large fixed /i. Might it be that the slow-roll conditions 
are satisfied with tan/3 — > 1? In fact it is easy to show that 

e „ (±?*LY „ ( 4 ^ /2 ) CQt2 ^ V ( 5 13) 

P \V D d{3 J \2X 2 sm 2 2f3+{g 2 + g' 2 )cos 2 2(3j K ' ' 

from which one sees that as (3 — > it /4, ~ [(g 2 + g' 2 ){jt — 4/3)/A 2 ] 2 . Thus we see 
that, unlike in the MSSM case, models with a sufficiently slow approach to the D-flat 
direction might be entertained. 
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6. Conclusions 



In conclusion, we have demonstrated that, in the presence of non- minimal gravitational 
interactions in the NMSSM, but not in the MSSM, there is contained in the Higgs sector 
a viable inflaton candidate, as in the SM. Although we have not yet completed an 
analysis of the running of the couplings in this model, we remarked earlier that the 
non-renormalization properties of SUSY and SUGRA models make the behavior of the 
non-minimal couplings different. In particular, unlike in the nonsupersymmetric case, 
the beta function for \ will therefore have a fixed point at % — 0. It would be interesting 
to implement refinements taking into account renormalisation group evolution described 



(and debated) for the SM case in Refs. [0-[|nj, and explore the consequences for the 



Higgs spectrum; in particular, perhaps in the MNMSSM |L9| where p = 0, the spectrum 
is more constrained. 

Given that these calculations seem to be concerned with the dynamics in the region 
beyond the GUT scale, one question is whether the same sort of results can be obtain 
within supersymmetric GUTS [p0|-[p3| such as SU(5) or S'O(IO). Does inflation end 
somewhere below the (reduced) Planck scale Mp ~ 2.4 x 10 18 GeV but above the GUT 
scale Mjj ~ 10 _2 Mp? Is the dynamics of inflation viable over such a relatively narrow 
range of energies as this? 

The most important unanswered question is whether these calculations make sense 
from the point of view of effective quantum field theory. As we discussed, this is a 
point of controversy in the literature, and we hope to return to this issue in the future. 
It is partly a question of the range of scales over which perturbation theory can be 
used, and it is partly a question of just what is the expansion parameter. We certainly 
have assumed that it makes sense to use the nonpolynomial potential on energy scales 
\f^h > Mp in the Einstein frame associated with changes in the gravitational constant 
in the Jordan frame, where the action has the familiar polynomial form. Even if the 
arguments were correct that, in the SM, this model for inflation is unnatural, will 
the situation be improved as a result of supersymmetry? We leave this for future 
investigation. 
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